Large N behavior of Non-Ohmic Tunnel Junctions 
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Non-Ohmic tunnel junctions are believed to occur in systems where the exciton and orthogonality 
catastrophe effects significantly modify tunneling rates. Here we present a simple non-perturbative 
treatment of the thermodynamic and transport properties of subOhmic and superOhmic tunnel 
junctions. Our analysis demonstrates the existence of a quantum phase transition in subOhmic but 
not superOhmic tunnel junctions. In addition, we find that the Coulomb gap vanishes continuously 
as the transition is approached. 
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Recently there have been several papers which have 
emphasized the importance of the Fermi-edge singular- 
ities to electron tunneling phenomenon. These include 
a treatment of Coulomb blockade oscillations in double 
quantum dots by Matveev and Glazman Q and the metal 
insulator transitions in granular materials by Drewes, 
Renn and Guinea M. In the later work, the authors 
argued that the metal-insulator transition observed by 
Herzog et al || in granular metallic wires was due to 
a quantum phase transition which occurs in sub-Ohmic 
tunnel junctions. Both these studies reflect the growing 
belief that exciton and orthogonality catastrophe effects 
are likely to be significant in a variety of systems in- 
cluding single tunnel junctions, quantum dot arrays, and 
granular metallic composite materials. 

The interest in these effects first began with the clas- 
sic treatments of x-ray absorption by Mahan Noziers 
and De Dominicus(MND) H,||. However, it was only 
a few years ago that the effect of Fermi edge singulari- 
ties on tunneling phenomenon was first considered. Per- 
haps the first work examine these effects was by Ueda 
and Guinea [@,|7| who showed that orthogonality catas- 
trophe effects could cause the conductance to behave as 
G ~ T x where x > 0. A subsequent study of reso- 
nant impurity tunneling by Matveev and Larkin B pre- 
dicted a zero temperature I(V) characteristic of the form 
I(V) oc (V - V t h)~ a 0(V - V th ). This prediction was later 
confirmed Q in an experiment involving resonant tunnel- 
ing between a 2-dimensional electron gas (2DEG) and a 
localized impurity. Based on these works it is now gener- 
ally believed that Fermi edge singularities a likely to pro- 
duce a variety of interesting effects including non-Ohmic 
behavior in the I(V) characteristics of tunnel junctions. 
For this reason, we will hereafter refer to such systems as 
non-Ohmic tunnel junction (NOTJ's). 

In order to better understand the behavior of NOTJ's, 
we will present a simple yet powerful non-perturbative 
treatment based on the large N expansion. Our approach 
provides convenient alternative to a renomalization group 
treatment of NOTJ's. In particular, one can easily cal- 
culate critical exponents characterizing the metal Insu- 
lator transition which occurs in some subOhmic tunnel 



junctions. In addition, one can also obtain detailed ex- 
pressions for conductances and specific heats. These cal- 
culations confirm the earlier prediction that sub-Ohmic 
tunnel junctions exhibit a quantum phase transition be- 
tween an insulating state to a state with a divergent zero 
temperature conductance. In addition, the calculations 
indicate that the insulating state exhibits a renormal- 
ized Coulomb pseudo-gap which collapses in a continu- 
ous fashion as the transition is approached. Finally, the 
large-N approximation indicates that the metallic state 
is absent in Ohmic tunnel junctions. 

In order to study the quantum dynamics of mesoscopic 
voltage fluctuations in NOTJ, will used a the long-range 
XY model introduced by Drewes et al. Q This model is 
similar to that first introduced by Ben Jacob, Mottola, 
and Schon ||To| to treat Ohmic tunnel junctions. The 
model is defined by the imaginary time partition function 



J Vn <5(n 2 - 1) exp -S[n] 



(1) 



where n(r) = (cos <^(r), sin(^(r)) may be used to relate 
the XY spin orientation to voltage using 4>{t) = eV/h. 
As discussed by Drewes et al 0], the effective action of 
NOTJ is 



S[n] = J drdr'air - t')[1 - n(r) • n(r')] 



(2) 



where ot(r) — cto T Q £ / r2_e ■ The parameter ao is a dimen- 
sionless parameter proportional to the squared tunneling 
matrix element. For e = 0, this model reduces to that of 
Ben Jacob et al. 

The above model has been analyzed Drewes et al using 
Monte Carlo simulations and the renormalization group. 
The renormalization group analysis shows that the model 
exhibits a correlation time 



£,t ~ \ao 



a c 



-l/e 



(3) 



where the phase transition occurs at a c = \/2ir 2 e. In the 
insulating phase one can associate h/£ T with a renormal- 
ized Coulomb gap subject to some important qualifica- 
tions: (1.) The large ao (spin- wave) calculation demon- 
strate dl/dV ~ T 2 ^^ behavior. So we expect that a 
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well defined Coulomb gap does not occur. However, at 
T = the leading order perturbation theory does give 
a threshold at Eq = e 2 /2C. We expect therefore that 
higher order cotunneling processes will put states within 
the gap. So A ~ Eq((x c — ao) _1 ^ e will n °t be a well 
defined gap although it may correspond to a pseudogap. 

The existence of pseudogap together with its behav- 
ior near the phase transition should be observable in the 
behavior of the non-linear I(V) characteristics. More 
generally, non-linear I(V) characteristics could provide a 
powerful means to test the scaling properties of the tran- 
sition. In particular, one expects such data according to 
the scaling ansatz scaling law 

dl/dV = & - (k B T/Ay> +1 -'F ± (eV/A, k B T/A) (4) 

where F_ (x, y) and F + (x, Y) are the a > a c and a < a c 
branches of a universal scaling function which is finite for 
x = y = 0. The V — > limit of this scaling form has been 
verified using Monte Carlo finite size scaling methods by 
Drewes et al al. 

A final consequence of the RG analysis is the relation 
rj = 1 + e ]l2| ]. This result is very important since it 
implies a temperature independent conductance at the 
subOhmic to insulator critical point. Although the con- 
ductance is obtained from the fixed point Hamiltonian 
and is unaffected by irrelevant operators, it is not uni- 
versal in any empirically useful sense. This unfortunate 
result follows from the fact that the model exhibits a line 
of fixed point Hamiltonians TL t which enable the "uni- 
versal conductances" to take any value ranging from to 
oo. (See discussion below.) 

At this point, we turn to our treatment of the large 
N approximation. We begin with a generalization of 
the long range XY model to a long range Heisenberg 
model. The dynamical degree of freedom of our model 
is now an N component unit-vector spin-field, h = 
(m, n-2, . . . , njv). Next we eliminate the n 2 = 1 con- 
straint by introducing an auxiliary field A(r) in eq. [l]. 
This gives 



Z = J T>[n,X] exp-J 



duj o:qCjv 
2tt 2 



iJdrX(r)[n 2 (r)-l} 



where C N = -4tqT(c - 1) sin(e7r/2). 

As N — > oo, the functional integral is dominated by 
a saddle point on the imaginary A(r) axis @. Small 
fluctuations about the saddle point will correspond to a 
theory with a finite correlation length £ T = h/A. Hence, 
it will be convenient, to rewrite the auxiliary field in the 
form A(t) = o:oCn{u(t) — iA 1 ~ e /2] Next we integrate out 
the n(r). This gives 



S effW) 



^(f -l)ln(CWao)- 



+ f In 



J draoCN(- 



— +m(r)) 
I- A 1_£ + 2iu(r)] 



We then obtain A by solving SS e f f[u]/5u(T)\ u=0 = or 
equivalently by solving 



This gives the 



where 
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A = E Q {\ - -i) 

a, 



[Eh~^] (5) 



(6) 



N 
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47re T(e - 1) sin(^f) 
The e — > limit, behavior of the A is 

27r 2 a 



A = E Q cxp-- 



N 



(7) 



(8) 



Next we consider spin-fluctuations about the saddle 
point. In the large N limit, G(t — r') =< n(r) • n(r') > 
is given by 



G(io) 



N 



1 



a Q C N \u\^ + A^ 



(9) 



This result demonstrates that £ T = h/A is indeed a 
correlation length. The fact that £ r diverges when 
do = a c then implies that a c — N/2ir 2 e is a critical 
point. This leading order in 1/JV result is, of course, 
consistent with the renormalization group results that 
a c = (N - l)/27r 2 e[l + 0(e)]. At the critical point the 
spin-spin correlation function becomes 



G(T)=r(e+l)cos(y)(Sy 



(10) 



This form implies i] = 1 + e as was previously obtained 
by Fisher Ma and Nickle jl2| and is believed to be exact. 

Next we wish to consider the behavior of the dimen- 
sionless XY model specific heat. The quantity under 
discussion is Cxy = a 2 l d 2 hiZ/da 2 ] . This is not the 
physical specific heat of a tunnel junction. It is, how- 
ever, a quantity which is relevant to Quantum Monte 
Carlo studies |1J] of the tunnel junction. The calcula- 
tion the heat capacity proceeds as follows. First, we cal- 
culate F = —ksTlnZ to leading order in 1/N using the 
u(t) = saddle point. This gives a zero temperature free 
energy 

F = fE Q HC N a ) - ^oCnA 1 ^ 

This coincides with kgTS[u = 0]. Now using the Eulcr- 
Lagrange equation for A together with eq. || one obtains 
the expression 
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2ir e 2 



(11) 
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This result is valid to leading order in 0(1/N) provided 
that ao < a c and \a c — ao\/ao are both small. One sees 
that the specific heat exponent a = 2 — 1/e coincides with 
that obtained from the Josephson relation a = 2 — dv. 

At this point, we wish to consider the behavior of the 
dc conductance. First we generalize the electric current. 
For the physical U(l) = 0(2) model, the current is ob- 
tained using J(t) = — 5S[A x }/5A x (t) where 



S[A X ] 



= /'drdr / a(T-r / )[l-cos($(T)-$(T'))] (12) 



where $(r) 



j-A x (t). This gives 



I(t) = 2e J dr'a(T - r') sin(^(r) - c/>(t')) 



(13) 



Equivalently, an expression for the electric current could 
be obtained by considering the variation of S under a 
local gauge transformation of the form <5</>(t) = 8e(r). 
Then ^8S/8c(t) — J(r). This second definition is con- 
venient for defining the electric current(s) in the O(N) 
model. In the O(N) model one has a set of N — 1 electric 
currents, I a (r) 1 a = 1 . . . N — 1 each associated with the 
N — 1 Lie algebra generators T a . Since n(r) is an N 
component vector field which transform according to the 
fundamental representation of O(N), the T a , matrices 
are N x N anti-symmetric matrices which can be cho- 
sen such that TrT a T b = 8 ab . We now define the electric 
currents as follows: Under a local gauge transformation 
8n(T)/8t a {r') = T a • u(t)5(t - t'). The currents are 
then defined to be I°(t) = f<5S[e]/ 'Se a (r). From this we 
obtain the explicit form of the electric currents 



I a (r) = 2e J dr'a(T - r')n(r) • T a • n(r') 



(14) 



The special case of the 0(2) = U(l) model, the above 
result reduces to eq. [l3| as required. This may be shown 
by taking T = ia y and n(r) = (cos(0(t)), sin(<^(r'))). 

Having defined the O(N) electric currents, we want to 
perform a Kubo calculation of the dc conductances G ab . 
The first step is to use the identity 



< 



ss ss 

Se a (r) 5e b {T') 



>- 



< 



5 2 S 



5e a (T)Se b (T') 



> 



(15) 



which follows from a functional integration by parts. This 
gives a current current correlation function for the phys- 
ical model (i.e. U(l) = 0(2) model) of the form p3| 

< I a (T)I b (r') >= -\a(T - t') < n(r) • T a T b • n(r') > 

(16) 

Using the Kubo formula and eqs. |l0| and |l6| we obtain 
the critical conductance, G c , of the tunnel junction to 
leading order in l/N: 



G ab = 2tt(1 - e)ctn 



-S a b = G c S a b 



(17) 



Because of the value of the r\ exponent the conductance 
is temperature independent. Next we observe that the 
conductance diverges as e — » 0. See fig. 1. This means 
that for small e a rather large value of ao (or equiva- 
lently the tunneling matrix element) is required in order 
to obtain a subohmic phase. This is consistent with the 
belief that the ordered (subOhmic) phase is destroyed as 
e — > and that the phase transition is absent for positive 
e. The second observation is that as e — > 1 a — > 0. This 
is consistent with the fact that the disordered (insulat- 
ing) phase is absent for e > 1. In particular, it implies 
that an arbitrarily weak ao will order the XY model as 
e — > 1 is approached. The fact that the critical conduc- 
tance can, depending on the value of e, range from to 
oo is interesting since it implies that no universal value 
of the conductance should be expected. 
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FIG. 1. The critical conductance vs. e. Observe that 
the insulating phase is squeezed out of the phase diagram 
(i.e. G c — > 0) as e — > 1. Similarly the conducting phase is 
squeezed out as e — > 0. 

Next we consider the AC conductance G(oj). For this 
purpose it is convenient to approximate < n(r) • n(0) > 
as being given by its critical behavior for r < < 1/ A and 
exhibiting a cutoff for r ~ 1/A i.e. we assume that 



G(r) «r(l + e )cos(y)^exp-(A|r|) 



(18) 



With this approximate form one can then obtain the 
time-ordered real time current current correlation func- 
tion: 



Cf(uj) = -t dt < T t (I a (T)I b (0)) > expiojt (19) 

J — oo 

Using eq. [l6| one obtains 
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Cf{w) = -2AG c S at 



duj' 



^ 2ir (to' - u>f 



16 



(20) 



This in conjunction with the the Kubo formula is then 
used to obtain the the ac conductance G(oj) ab = [G'(u>) + 
iG" (ui)]6 a b- The real part is found to be 



G'H = G c (l 



2\oj 



■)0(\uj\ - A) 



(21) 



where 6{x) — for x < and 1 for x > and where G c is 
the critical conductance as given by eq. [l?]. This result 
allows one to identify A as an excitation gap which, near 
the transition, is much smaller than the charging energy 
Eq- 

Previously, Drewes et al H found a well defined 
Coulomb gap in a leading order ( O(ao)) perturbative 
calculation of the /(V) characteristic. However a a 2 , 
culation of G(T) found that 



cal- 



G 



al ( nk B T 



2(l-£) 



(22) 



That result indicates that the conductance does not 
exhibit an activated temperature dependence and that 
higher order processes introduce states inside the gap. 
Moreover, one would argue that the zero temper- 
ature I(V) characteristic is of the form I(V) ~ 
(e 2 /h)al(eV/A) 2 ^^ for V « A and that ac conduc- 
tance is of the form 



G'(u) ~ (w/A) 2 ^ 



(23) 



when uj « A. These results presumable occur at higher 
than leading order in the 1/N expansion. 

In addition to calculating the G'(uo) one can also calcu- 
late the imaginary part of the conductance. In this case 
one finds 



Ct(0) 



—G c 

TTOJ 



In 



A 2 - . 



A ln 



(24) 



where Ct(0) depends logarithmically on Eq/A. 
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